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Coupled quantum electrodynamics (QED) cavities have been recently proposed as new systems 
to simulate a variety of equilibrium and non-equilibrium many-body phenomena. We present a brief 
review of their main properties together with a survey of the last developments of the field and some 
perspectives concerning their experimental realizations and possible new theoretical directions. 
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I. INTRODUCTION 

The ability to design and fabricate controllable many- 
body systems has been realized to be a precious tool 
to explore the world of strongly correlated systems. In 
several different physical phenomena, ranging from high 
temperature superconductivity to heavy fermions or the 
fractional quantum Hall effect, strong local electronic 
correlations play a crucial role. The seemingly simpli- 
fied models to describe these correlations, the Hubbard 
model of high temperature superconductivity for exam- 
ple, are extremely difficult to solve and over the years 
a number of judicious analytical and numerical methods 
have been developed. Together with these more "tra- 
ditional" methods, in recent years it was proposed that 
strongly correlated systems could be studied by means of 
quantum simulators [T], i.e. fabricated systems that can 
experimentally simulate the model Hamiltonian underly- 
ing the non-trivial properties of the physical systems un- 
der consideration. The advantages of this approach are 
twofold. First of all it is possible to explore the properties 
of strongly correlated model Hamiltonians also in those 
regions of the phase diagram which are elusive to numer- 
ical and analytical investigations. Secondly it allows to 
test to which extent the model Hamiltonians under con- 
sideration are appropriate to treat the physical systems 
that they are supposed to describe or whether additional 
ingredients are necessary. 

Quantum simulators have a long and successful story. 
Probably the first fabricated system to have these char- 
acteristics were Josephson junction arrays Q\. The field 
boosted with the appearance of cold atoms in optical lat- 
tices [3] which proved to be excellent simulators of a large 
variety of strongly interacting Fermi and Bose systems. 
Topic of the present paper is to review the latest develop- 
ments of a newly born direction in the field of quantum 
simulators based on arrays of QED-cavities [IHSj • There 
is already a comprehensive review on the subject [7], in 
this paper we will focus on some of the recent develop- 
ments. 

Cavity-QED arrays offer the possibility to realize 
strongly correlated states of light. They can operate at 
high temperatures (as compared to Josephson arrays and 
optical lattices) and allow for single-site addressing thus 
opening a way to access experimentally correlation func- 
tions. Furthermore they might allow to explore a number 
of new equilibrium and non-equilibrium quantum phase 



transitions. 

This paper is organized as follows. In the next section 
we will introduce cavity arrays. We define the model 
Hamiltonians that govern their dynamics and the most 
relevant sources of dissipation. In Sec.|Tn]we discuss the 
most important characteristics of the equilibrium phase 
diagram and the nature of the low-lying excitations. Sec- 
tion [IV] describes how cavity arrays can operate as quan- 
tum simulators. Section [V] reviews some properties of 
these systems out of equilibrium. We complete the pre- 
sentation with the conclusions where we will briefly dis- 
cuss possible experimental realizations and some perspec- 
tives in the field. 



II. THE MODELS 

A sketch of a cavity array is illustrated in Fig. [T] It 
consists of a regular arrangement of QED-cavities which 
are coupled by the hopping of photons. Light resonates in 
each cavity and interacts with matter contained therein. 
Experimental realization of cavity arrays can be imagined 
in different ways. We will briefly discuss this point at the 
end of the paper. Fig. [l] illustrates the case in which the 
array is realized using photonic crystals. The salient in- 
gredient at the base of the rich physics of cavity arrays 
is the interplay of two competing effects. Light-matter 
interaction inside the cavity leads to a (possibly strong) 
nonlinearity between photons. On the other hand, pho- 
ton hopping between neighboring cavities favors delocal- 
ization thus competing with photon blockade. A model 
that describes a cavity array must then take into consid- 
eration the interaction of the light and the matter within 
each cavity, the coherent coupling between the cavities 
induced by the propagation of the light, the leakage of 
photons out of the cavities, and the dissipation and de- 
coherence of the matter. It is also necessary to include 
an external pump to populate the cavities and devise 
a method to perform the measurement of the state of 
each cavity. We first describe the terms that lead to the 
unitary evolution of the array and afterwards the main 
sources of dissipation. 

A single cavity confines several modes of the electro- 
magnetic field and each mode is quantized as an harmonic 
oscillator. In the case that a single mode with frequency 
wc is considered, the corresponding Hamiltonian is given 
by He = i^a^iat where the operator an (a|) annihilates 
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The Jaynes-Cummings model reads 



hopping 



FIG. 1: A sketch of a QED-cavity array. It consists of a 
regular arrangement of QED cavities. Neighboring cavities 
are coupled by photon hopping. Nonlinearities in the cavi- 
ties may produce an effective repulsion between the photons 
leading to an anharmonic spectrum. The nonlinearity may be 
produced, e.g. by a two-level system (depicted in the inset) 
coupled to the light resonating in the cavity and subjected to 
decay. Photons in the cavities have a finite lifetime therefore 
the cavities are pumped with an external coherent drive. 



(creates) a quantum of light in the mode of the ^-th cav- 
ity. If the cavities are sufficiently close to allow for photon 
hopping, an additional kinetic term — J(aJ,a£-|-H.c.) (J is 
associated to the tunneling rate) should be added to the 
Hamiltonian. These two terms constitute the model for 
a cavity array in the absence of any interaction with the 
matter field. In the presence of hopping the Hamiltonian 
of the photons is still harmonic and can be diagonalized 
by Fourier transform leading to the dispersion law of the 
photons in the lattice. 

The presence of atoms (quantum dots or Cooper pair 
boxes depending on the implementation) inside each cav- 
ity leads to a strong effective nonlinearity between pho- 
tons. It is enough to think to the matter field as a few- 
level system coupled to cavity mode and possibly to some 
external source. The Hamiltonian for an array of cavities 
can be written as on each cavity and the photon hopping 
term between different cavities 



^ = H '^^"^ ~ -^ H (^^"^' + H-*^-) 
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The local contribution T-L^ ' describes the light-matter 
interaction. In the limit in which leakage of photons is 
ignored, the Hamiltonian in Eq. M has been treated in 
the grand-canonical ensemble. In this case one should 
add a term containing the chemical potential. In the 
next section the phase diagram will be presented in the 
grand-canonical ensemble. 

Probably the simplest model to describe the interac- 
tion between light and matter is the Jaynes-Cummings 
model in which one mode of the cavity interacts with a 
two-level system [8j . One photon can be absorbed by the 
two-level system, that goes into the excited state, and 
conversely a photon can be emitted if the two-level sys- 
tem flips from the excited |2) to the ground state |1). 
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a g are the raising/lowering operators for the two-level 
system and e denotes the transition energy between the 
two levels. In the rotating frame with respect to the un- 
coupled Hamiltonian the relevant quantity is the detun- 
ing A = w — e. The spectrum of Eq. (pi) is anharmonic so 
that, effectively, the two-level system induces a repulsion 
between the photons in the cavity. The strong effective 
nonlinearity between the photons turns the cavity into 
a turnstile device, where only a photon can be present 
at the same time. Intuitively, this can be understood as 
the fact that one photon in the cavity strongly modifies 
the effective resonance frequency, inhibiting the injection 
of a second photon. This phenomenon has been termed 
"photon blockade" 15^ , after the Coulomb blockade effect 
of electrons in mesoscopic structures. The entire cav- 
ity thus behaves as an effective spin system, that emits 
strongly anti-bunched light. The cavity array Hamilto- 
nian with the Jaynes-Cummings term was first used by 
Greentree et al. in Ref. [S] and by Angelakis et al. in 
Ref. [5] . A similar model with many two-level systems in 
each cavity has also been used [TMT^ . where each two- 
level system interacts with the mode of the cavity in a 
sort of Dicke-Bose-Hubbard model. Na et al. in Ref. [TT] 
allowed for the number of systems to fluctuate among 
the cavities and showed that if the cavity is strongly red- 
detuned with respect to the transition in the two-level 
systems, the photons experience an effective Kerr non- 
linearity. Koch and Le Hur |13) . though, pointed out 
that the substitution of the Jaynes-Cummings coupling 
with an effective Kerr nonlinearity is not appropriate in 
general. 

Since the quantum effects in the cavity may be visi- 
ble if the coupling exceeds the decay amplitude of the 
matter and of the light, there has been a lot of effort 
in increasing the effective cavity nonlinearity. Strong 
interaction of light is in general limited by the absorp- 
tion of the medium, because the optical nonlinearities 
are weak in non-resonant processes, while the absorp- 
tion plays a dominant role in the resonant processes. 
One route to access the strong-coupling regime involves 
the usage of external classical sources that prepare co- 
herently the matter in a state with reduced absorption. 
Electromagnetically Induced Transparency (EIT) can be 
used to achieve interaction strengths that are one order 
of magnitude larger than what is possible with an ensem- 
ble of two- level systems [T3] . Various schemes to obtain 
strong nonlinearities have been discussed [HI fT5HT7] . In 
this spirit Hartmann et al. in Ref. [4] considered an en- 
semble of four-level atoms in a cavity array defined, in 
the rotating frame, by a local interaction of the form 
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+ 31 (5f a\ + Sf at) + 52 (5f al + Sf at) , 
having defined the global atomic raising and lowering 



operators 5'™ = Sj=i IOj("^li f'''™ — IOi(™l ^re the 
atomic raising and lowering operators {I ^ m), or en- 
ergy level populations (/ = m) for the jth atom. The 
transition |3)j — ?> |2)j is driven by a classical coupling 
field with Rabi frequency J7; the cavity mode couples the 
\l)j -^ \3)j and \2)j — )■ |4)j transitions with coupling con- 
stants gi and 52; the parameters (5 and A account for the 
detunings of levels 3 and 4, respectively. They considered 
a setup similar to Ref. [9] and introduced three families of 
polaritons, delocalized over the whole atomic ensemble, 
that diagonalize the Hamiltonian (if the state |4) decou- 
ples from the dynamics and the number N of atoms is 
very large). One of the polaritons projects only on the 
lowest- lying metastable states |1) and |2) and is thus the 
longest-lived. The polaritons can be considered indepen- 
dent bosonic particles if their number is much smaller 
than the number of atoms inside each cavity [15]. The 
longest-lived polariton decouples from the others and ex- 
periences an energy shift due to the perturbative coupling 
to the level |4), that induces an interaction of the Kerr 
form. In this way, in the large N limit, they were able 
to show that the cavity model maps onto an effective 
Bose-Hubbard model [11] for polaritons 
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Up to now we did not take into account the leakage 
of photons out of the cavities and the decoherence and 
dissipation of the few-level systems. One possibility, in 
the case of slow decay, is to complement the Hamilto- 
nian with imaginary frequencies multiplied by the pro- 
jectors of the decaying states, that render the dynamics 
non-Hermitian (see some application of this method in 
Ref. [14 ). Alternatively the decay can be given in terms 
of the density matrix p of the cavity array, using the 
master equation 

dtp{t) = -i[H,pit)]+C[pit)], (5) 

where the Liouvillian C in the Lindblad form reads e.g. 

^Pit)] ^ 2 ^ {at,p{t)al - alaip{t) - p(t)a|af j . (6) 

e 

This term describes leakage of photons from all the cav- 
ities of the array, with an equal rate k. A similar term 
can be used to describe decoherence and dissipation of 
the matter, although these processes can usually be ne- 
glected in the timescale of the photonic decay [7] . 

In order to have a measurable signal out of the cavity 
array it is necessary to refill with photons the modes of 
the cavities, to contrast the leakage. If a coherent laser 
beam of frequency wl is coupled to the cavity, the Hamil- 
tonian acquires a term e{t)Lal + e{t)*j^ai, where e(i)L is 
proportional to the electric field of the beam. To elimi- 
nate the explicit time-dependence e*"^* from e(i)i, it is 
convenient to use a rotating frame in which the strength 



of the pumping is constant and the frequency wl is sub- 
tracted from the energy of the photons in the cavity 
mode. The contribution of the drive to the Hamiltonian 
reads 



ni.B 
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In the next section we will discuss the main features of 
the phase diagram deriving from the Hamiltonian (flj). 



III. PHASE DIAGRAM 

The phase diagram of the Bose-Hubbard model defined 
in Eq. Q has been studied for more than two decades 
(see for example references in [3j). A qualitative under- 
standing of the zero-temperature phase diagram can be 
obtained by considering the two limiting cases in which 
one of the two coupling energies (J or U) dominates. If 
the hopping dominates, the bosons delocalize through- 
out the whole lattice. In the limit of vanishing inter- 
actions, the many-body ground state becomes simply an 
ideal Bose-Einstein condensate where all the bosonic par- 
ticles are in the Bloch state of the lattice with vanishing 
quasimomentum, so that the many-body wavefunction is 



written as I -7= X]f=i 'A 



N 



Ivac) , where L is the total 



number of cavities, N is the total number of bosons and 
|vac) is the vacuum state of the many-body Fock basis. 
In contrast, if the interaction dominates, each site has a 
well defined number of bosons in the ground state. In 
order to put an extra boson on a given site, one has to 
overcome an energy gap of the order of U. In the Mott 
lobes each cavity contains a definite number of bosons, 
the fluctuations of the number operator vanish, and each 
particle is localized in a given cavity. Each Mott lobe 
exists up to a maximum hopping (J/C/)crit at which it 
is energetically convenient for the particles to delocalize 
over the whole lattice (decreasing their kinetic energy) 
rather than localize in a given site (decreasing their mu- 
tual interaction). In the case of vanishing hopping, the 
Hamiltonian is the sum of commuting local Hamiltoni- 
ans and the ground state can be written in the product 
form Yli \n)e- When the hopping is comparable with the 
onsite repulsion a quantum phase transition takes place, 
separating a superfluid from a Mott insulator. Although 
there are large quantitative differences between the case 
of one-, two-, or three-dimensional lattices, the equilib- 
rium phase diagram of the model (at zero temperature, 
in the plane of the hopping and the chemical potential) is 
similar in all cases and consists of a superfluid phase sur- 
rounding a sequence of Mott-insulator lobes [see Fig. [2] 
(left panel)]. The U{1) symmetry is broken in the super- 
fluid phase, where the kinetic energy dominates. 

The phase diagram of the cavity array described by 
Eq. (IT]) resembles strong similarities with that of the 
Bose-Hubbard model. In Fig. [2] (right panel) we report 
for comparison the phase diagram obtained in Ref. |10| 
for a one-dimensional array. Quantitative differences 
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FIG. 2: Left: sketch of the phase diagram of the BH model 
Q. The green regions are the Mott lobes, z is the coor- 
dination number of the lattice. Each lobe extends up to a 
maximum critical ratio {J/U)ciit- Right: the phase diagram 
of the cavity array in the case of a Jaynes-Cummings model, 
according to the DMRG computation presented in Ref. [10] . 
The green area corresponds to a Mott insulator of polaritons 
with average filling n — 1, 2, and 3. 



arise depending on the model of the cavity. For the 
sake of clarity we will ignore these differences and try 
to highlight only those characteristics which are common 
to cavity arrays. The phase transition between the Mott- 
insulator and the superfluid is a symmetry breaking phe- 
nomenon that involves the conservation of the number 
of bosons J2e''^i- The conserved quantity in the cavity 
model, take for example the Jaynes-Cummings interac- 
tion, instead is the number of polaritons, that do not 
have a definite quantum statistics. In cavity arrays the 
detuning A can also be varied so that the phase diagram 
becomes effectively three-dimensional [T^. It turns out 
that the detuning is a convenient experimental knob to 
tunc the system across the phase transition. Moreover, 
most importantly, to change the detuning from negative 
to positive values drives the nature of the excitations 
from quasilocalized excitons to polaritons to weakly in- 
teracting photons [20 . 

A cavity array exhibits spectral properties similar to 
the BH models, including gapped particle and hole bands 
in the Mott insulating phase and Bogoliubov-type exci- 
tations in the superfluid phase. The single-particle exci- 
tation spectrum in the Mott phase has been investigated 
also by Aichorn et al. |20j using a variational cluster ap- 
proach, in which the self-energy of the infinite lattice 
of cavities is approximated by a finite reference system. 
Schmidt and Blatter |2T], for the Jaynes-Cummings case, 
found four modes corresponding to particle or hole exci- 
tations, versus the two modes of the BH model. Pippan 
et al. in Ref. |22] studied the dynamical structure factor 
and the single-particle spectrum using a quantum Monte 
Carlo simulation in one dimension. An analysis of the 
nature of the excitations in small clusters was also per- 
formed in Ref. [53]. Paternostro et al. in Ref. [23] inves- 
tigated the dynamical properties of an array with many 
two-level systems in each cavity, beyond the elementary 
excitation picture and found that the time-evolution of 
the array, in the mean-field approximation, supports a 
soliton-like excitation. 

Following the initial investigations [IHS], the existence 



of the Mott lobes in the cavity arrays has been verified 
with several different methods. Makin et al. [25] resorted 
to the exact diagonalization of small sets of cavities with 
several topologies and to a cluster-mean-field approxi- 
mation. Schmidt et al. in Ref. [21 proposed a strong- 
coupling theory based on a linked-cluster expansion of 
the Green's functions at finite temperature. This method 
includes quantum fluctuations beyond the mean-field re- 
sults, that are recovered instead in the random-phase ap- 
proximation. With these spectral methods the border of 
the Mott lobe is signaled by the vanishing of the gap. 
A quantum Monte Carlo analysis in two dimensions has 
been performed by Zhao et al. in Ref. |26j . The quan- 
tum Monte Carlo results for the phase border compare 
quantitatively with the analytical solution [21 almost up 
to the tip of the Mott lobe. In one dimension, the exis- 
tence of the insulating Mott lobes has been established 
in Ref. [lO] by means of a DMRG calculation. With the 
same method, it has been shown in Ref. [57] that the 
polariton fiuctuations in a finite block of cavities yields a 
clear signature of the two phases. The universality class 
of the transition appears of the same type as the BH 

[111121]. 

The nature of the Mott-insulator for polaritons is dif- 
ferent from the BH case, in that the number of polari- 
tons is fixed to an integer value, but the fluctuations of 
the number of photons (ne) do not vanish. Similarly, 
in Ref. [22] it was shown that the exciton and photon 
structure factors cannot be used to characterize the Mott 
state. Although the photons are not the genuine excita- 
tion of the system, it is desirable to characterize the tran- 
sition in terms of photonic observables since the emission 
of light from the cavity can be measured in experiments. 
To measure the polaritonic state using the projection on 
the atomic component [31 [6] it is necessary to address the 
system with further optical probes. 

The presence of many two-level systems in each cav- 
ity and its effect on the phase diagram has been ana- 
lyzed in [ini dl] , conflrming that even with a moderate 
number of atoms in each cavity there is a good quantita- 
tive agreement with the mapping onto the BH performed 
in [3] . New features appear if the number of atoms fluc- 
tuates from cavity to cavity. In this case a glassy phase 
for polaritons is expected around the Mott lobes |10| . 
Also disorder in the distribution of the excitation and 
cavity frequencies [llj , and the effects of a finite temper- 
ature [ini Hi] have been considered. 

Using more complex models of cavity arrays a two- 
component Bose-Hubbard model [55] and superradiant 
Mott phase [5S] have been shown to be realizable. 



IV. CAVITY ARRAYS AS QUANTUM 
SIMULATORS 



Building upon the models presented in Sec. |TT| in spe- 
cific parameter ranges or with the addition of further 
optical tools, it is possible to use the cavity array as an 



implementation of other many-body models. In other 
words, the cavity array offers an implementation of a 
"quantum simulator" for lattice models, in the same 
spirit as the Josephson junctions arrays or the cold atoms 
in optical lattices. The first requirement for a quantum 
simulator is to act as a calculator specifically tailored 
to the solution of the model that it implements. The 
knowledge on the model gained with the quantum sim- 
ulator can then be applied to all the other physical sys- 
tems described by the the same model. With respect to 
the implementation of a quantum simulator in an opti- 
cal lattice, the cavity array may offer the advantage that 
each site of the array can be addressed independently. 
The expectation values of the observables could then be 
measured directly from the light emitted by one cavity. 
Moreover, the cavity arrays can be implemented with 
several different experimental systems and this may offer 
some advantages. On the other side, optical lattices seem 
unbeatable in terms of scalability and absence of imper- 
fections. It should be said that at present there are yet 
no experimental realizations of cavity arrays. The hypo- 
thetical advantages of these quantum simulators need to 
be tested against experimental realizations. 

The earliest family of lattice many-body systems that 
have been investigated in the context of cavity arrays is 
spin lattice systems. The indication that a single cavity 
could act as an effective spin system for the photons dates 
back to the recognition of the photon blockade effect [S] . 
Angelakis et al. [6] derived an effective XY Hamiltonian 
from the Jaynes-Cummings model after decoupling the 
upper and the lower polariton in the Mott polaritonic 
phase. The realization of a ZZ Ising coupling ulcrl,-^ has 
also been discussed f30^ and the possibility to achieve 
both XY and ZZ couplings has been put forward by Hart- 
mann et al. in Ref. [31] . The proposal relies on a cavity 
array containing three-level atoms, in a A configuration. 
The two lowest-lying states |1) and |2) are coupled to 
the excited state |3) both via the quantized mode of the 
cavity and via external classical sources. With an ap- 
propriate choice of the couplings and the detunings, the 
dominant Raman transition between the two lowest-lying 
states involves one laser and one cavity photon. This 
implies that the emission and absorption of virtual pho- 
tons in the cavity is accompanied by a transition between 
the two lowest-lying states, that become the effective 1/2 
spin in the cavity. The coupling of the virtual photons 
between neighboring cavities produces a XY model. This 
effective spin can be coupled to an effective magnetic field 
in the z direction. With the same atomic spectrum but a 
different configuration of the external sources, it is pos- 
sible to implement a ZZ interaction of the form cr|CT|, j^. 
The XY and the ZZ couplings do not act at the same 
time, but are combined in a unique effective Hamiltonian 
with the application of the Suzuki- Trotter decomposition 
in which the lasers that produce either interaction are 
periodically flashed one after the other. Cho et al. [32] 
focused on the implementation of high-spin Heisenberg 
models. Given the lack of analytical or numerical infor- 



mations that are available today on the phase diagram of 
high-spin systems, the possibility to measure the phase 
diagram directly from a quantum simulator is very in- 
teresting. The scheme includes an ensemble of three- 
level atoms in a configuration similar to Ref. , 31j . but 
relies on fixed external lasers. The model allows to simu- 
late terms in the Hamiltonian of the form Sj, (5'|)^, ^f, 
'-'^ *-'^-j_i ! "-^^ '-'^-i-i ! '-'£*-'^j_i5 and o^o^,-^. l^uite a lew ex- 
ternal sources are necessary to control independently all 
the constants in the model, while the magnitude of the 
spin is determined by the number of two-level systems in 
the cavity. 

Taking into account the polarization of light in the mi- 
crocavities, Ji et al. in Ref. |33j put forward a model that 
is able to exhibit ferromagnetism of the photonic circular 
polarization. The model considers an hexagonal photonic 
crystal with a square superlattice of band gap cavities, 
doped with three-level systems. In the same system, it 
is also possible to access adiabatically the antiferromag- 
netic phase, and an exotic super-counter-fluid phase in 
which the total current of photons is zero as in the Mott 
phase, but the currents with deflnite circular polarization 
are nonzero and opposite. 

The local addressability of cavity arrays has been ex- 
ploited by Cho et al. in Ref. j^ to propose the simula- 
tion of hard-core bosons with Abelian vector potentials. 
In particular this method allows the simulation of the 
Laughlin wave function for the fractional quantum Hall 
effect. The proposal assumes a two-dimensional array 
of cavities, each containing a single three-level system. 
To address independently the two directions of the array, 
two modes of each resonator are used, with the assump- 
tion that the frequency difference between the two modes 
is much larger than the coupling strength g. The simula- 
tion of the hard-core bosons follows from the simulation 
of a spin system using the lowest-lying states of the A 
configuration, along the lines of Ref. [32]. The ampli- 
tude of the interaction between the spins on neighboring 
sites is proportional to a position-dependent phase, that 
defines the gauge potential. This phase is controlled by 
the external laser sources acting on the three-level system 
and the local addressability of the cavities is invaluable 
to the precise implementation of this phase. Moreover, 
by an adiabatic change of the laser phases, it is possi- 
ble to insert a flux quantum through the two-dimensional 
plane, thus creating and moving a quasi-excitation in the 
system. 

Cavity arrays can be used also to store and manipu- 
late the resources necessary to perform quantum com- 
putation. The flrst study in this direction has been 
made by Angelakis et al. in Ref. |35j with the aim to 
realize quantum gates with photonic-crystal waveguides. 
In this case the coupling between the photons and the 
matter is not present in each cavity, but is used only in 
some selected points of the quantum circuit to induce 
a phase gate between the qubits. The coherent control 
of the photonic transmission through a cavity array has 
been considered by Hu et al. in Ref. [36', where it was 




FIG. 3: Non-equilibrium signature of the quantum phase 
transition between the Mott-insulator and the superfiuid state 
in the Bose-Hubbard model Q with the leakage of photons 
described by (Im, obtained by Tomadin et al. in Ref. [40j. The 
order parameter ip = (a), averaged in an interval of time, is 
shown for [//«; = 50 (empty circles), 100 (filled circles), and 
200 (empty triangles) . The transition of the order parameter 
takes place together with a transition from antibunched to 
bunched light. On decreasing the dissipation strength the 
onset at which the order parameter becomes non-zero ap- 
proaches the equilibrium value (indicated as a vertical line). 



shown, by a Green's function technique, that the disper- 
sion of the hybrid hght-matter excitations in the array 
can be controlled by acting on the population inversion of 
the dopants. A procedure to generate indistinguishable 
single photons or polarization-entangled photonic pairs, 
which are important resources to the implementation of 
the quantum computation algorithms, has been demon- 
strated in Ref. |37j by Na et al. resorting to the exact in- 
tegration of a small number of cavities. The complete im- 
plementation of the Grover's search algorithm has been 
proposed recently by Kyoseva et al. |38l together with 
a proposal to implement one-way quantum computation 
with cavity arrays |39) . 



photons in a circular array of cavities in the Tonks- 
Girardeau limit. All the cavities are driven by a single 
laser beam and the far-field emission, that gives access 
to the occupation of the single-particle states in the mo- 
mentum representation, is computed. The laser couples 
to the many-body wavefunctions supported by the array 
depending on the total energy of the state. This yields 
a spectroscopic analysis of a many-body system that is 
created by the probe beam itself. In the limit of very 
strong Kerr nonlinearity, the photons become hard-core 
bosons and, being in one dimension, are represented by 
an effective wave function for noninteracting fermions. 
The classification of the states given by this mapping is 
used also to characterize the absorption in the case of 
intermediate interaction strength. 

The quantum phase transition of the BH (see Sec. Ill I, 
in an open cavity array, has been investigated by 
Tomadin et al. in Ref. "40] who proposed to discriminate 
the Mott-insulator and the superfiuid phases analyzing 
the light emitted by the cavity array following a pulse 
that creates a Mott state in the system. The existence of 
the two phases is clearly seen in the coherence properties 
of the emitted light (see Fig. Isl) for currently achievable 
values of the ratio between the interaction and the dissi- 
pation. The method relies on the dynamical instability 
of the BH that follows a quantum quench from the Mott 
to the superfiuid phase. The realization of the quantum 
quench in the dissipative environment is realized with a 
careful design of the initial pulse, that is substantially 
different from a 7r-pulse, and it is shown that the detec- 
tion of the phase border is very robust to imperfections in 
the pumping. The time-evolution of the system produces 
antibunched light if the ratio J /U is smaller than the 
critical value ( J/C/)crit at the tip of the Mott lobe, and 
the crossover from antibunched to bunched light clearly 
marks the phase boundary. 



V. NON-EQUILIBRIUM BEHAVIOR 

In order to detect experimentally the different quan- 
tum phases presented in Sec. |III| one assumes that the 
mixed light-matter excitations have a lifetime that is 
longer than the time necessary to the measurement pro- 
cess. Similarly, in order to perform the quantum simula- 
tions of a ground state discussed in Sec. |IV| it is necessary 
that the timescale of the many-body dynamics is shorter 
than the decay time. For example, the ground state of a 
many-body system could be reached if the decay of the 
polaritons is at least slower than the intercavity tunneling 
[TT] . As the decay times of the photons and of the matter 
strongly depend on the specific implementation [7], it is 
interesting to consider schemes that allow to measure the 
signatures of the many-body physics without relying on 
a negligible dissipation. It is then necessary to study the 
cavity array under strong non-equilibrium conditions. 

Carusotto et al. in Ref. [41, investigated a system of 



VI. CONCLUSIONS 

As candidates for simulating strongly interacting mod- 
els, coupled cavities present new characteristics as com- 
pared to other successful examples, like optical lattices 
or Josephson junction arrays: most notably it is possible 
to access their local properties. In this paper we briefly 
review their basic properties and some of the latest de- 
velopments in the field. 

QED-cavities can be realized in a number of different 
ways [1^44) . This flexibility in the design is a potential 
advantage to realize different local nonlinearities. More- 
over different systems may allow for different measure- 
ment schemes as well. As already mentioned, at present 
there are no experimental realizations of cavity arrays 
although it seems that technological requirements are al- 
ready at hand to fabricate a small number of coupled 
cavities. In this respect it may be interesting to explore 
properties of small clusters. 

In this review we conflned the discussion to coupled 



cavities. There are other systems in which many-body 
physics with hght can be reahzed (see for example [45l 
35]). Due to space hmitations we did not touch these 
very interesting directions of research. 
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